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Abstract 

We illustrate the use of the computer algebra system Macaulay2 
for simplifications of classical unirationality proofs. We explicitly treat 
the moduli spaces of curves of genus g=10, 12 and 14. 

If a moduli space M is unirational, then ultimately we would like to ex- 
hibit explicitly a dominating family, which depends only on free parameters. 
Although this is possible along the lines of a known unirationality proof in 
principle, this is far beyond what computer algebra can do today in most 
cases. 

However, if we replace each generic free parameter with a random choice 
of an element in a ground field F, then in many cases the computation of a 
random element of the family defined over F is possible. 

This is particular interesting over a finite field F, because in that case 
there is no growth of the coefficients in Grobner basis computations. The 
equidistributed probability distribution on F induces a probability measure 
on M(F) via the algorithm, and opens up for us the possibility to investigate 
the moduli space experimentally. With semi-continuity arguments we can 
sometimes deduce effectiveness of certain divisors on the moduli space. The 
computation with a single random example over a finite field can verify that 
the constructed family indeed dominates the desired moduli space. The key 
advantage in using Computer algebra instead of theoretical arguments is that 
random choices give smooth varieties almost certainly, while smoothness is 
always a difficult issue in any theoretical treatment. 

In this note we illustrate this technique in a number of cases, and provide 
explicit code for the computer algebra system Macaulay2. A more elaborate 
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code is provided with the Macaulay2 package RandomCurves, which will be 
freely available from |M2] and |BGS] . I recommend running the code available 
from www.math.uni-sb.de/ag/schreyer/home/forHandbook.m2 parallel with 
reading the code in the article. My favorite set-up of Macaulay2 is within 
emacs, which for example has syntax highlighting. That makes the code 
much more readable. Instructions how to use Macaulay2 with emacs come 
with the installation package of Macaulay2, see |M2j . 

Acknowledgements. I thank the referees for their careful suggestions, and 
Florian GeiB and Hans- Christian Graf von Bothmer for valuable discussions. 

1 Random Plane Curves 

Let DJlg denote the moduli space of curves of genus g. For a general smooth 
projective curve C of genus g the Brill-Noether locus 

W r d {C) = {Le Pic d (C) | h°(L) > r + 1} 

is non-empty and smooth away from W d +1 (C) of dimension p if and only if 
the Brill-Noether number 

P = p(g, d, r) = g - (r + l)(g - d + r) > 0, 

equivalent ly, iff 

d>{{r + l){g-r)-g)/{r+l). 

Moreover, W£(C) is connected in case p > 0, see |ACGH] . The tangent space 
at the point L e W£(C) \ W r d +1 {C) is naturally dual to the cokernel of the 
Petri map 

H°(C, L) ® H°(C, io c ® L- 1 ) -> H°(C, u c ). 

To prove the unirationality of S0T 9 we might prove that the Hilbert scheme 
of the corresponding models of C in P r is unirational. 

Severi's proof of the unirationality of S0T 9 for g < 10 is based on the fact 
that a general curve up to genus 10 has a plane model with double points 
in general position, see |Sev21[ IAS79] . An easy dimension count shows that 
this cannot be true for g > 11. 

The algorithm that computes a random curve of genus (7 = 10 proceeds 
in four steps. 

1. Compute the minimal degree d such that p(g, 2, d) > 0. 
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2. Choose a scheme A of 5 = ( d 2 1 ) — g distinct points in P 2 . 

3. Choose a random element in / e H c '(P 2 ,T\(d)) . 

4. Certify that / defines an absolutely irreducible 5-nodal curve. 

In case g = 10 we have d — 9 and 5 = 18. 

The Macaulay2 code for these steps looks at follows. The code is, I 
believe, fairly self explaining, since the Macaulay2 language is close to stan- 
dard mathematical notation. For further explanations I refer to the online 
documentation http:/ /www. math. uiuc.edu/Macaulay2/. 

Step 1. 

11 : — function that computes the minimal d: 

dmin= (r , g) ->ceiling( ( (r+1) * (g+r) -g) / (r+1) ) 

12 : g=10,r=2,d=dmin(r,g) ,delta=binomial(d-l,2)-g 

13 : rho=g-(r+l)*(g+r-d) 
o3 = 1 

Step 2. We specify 5 points in P 2 : 

14 : FF=ZZ/10007 — a finite ground field 

15 : S=FF[x_0. .x_2] — the homogeneous coordinate ring of P2 

16 : — pick a list of delta ideals of random points: 

point s=apply (delta, i->ideal (random (1 , S) , random (1 ,S) ) ) 

Remark: The Macaulay2 function "random(n,S) " picks a random element 
of degree n in a graded ring S. 

17 : IDelta=intersect points 

18 : betti res IDelta 

12 

08 = total: 14 3 

0: 1 . . 

1: . . . 

2: . . . 

3: . . . 

4: . 3 . 

5: . 1 3 
08 = BettiTally 
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For a very small finite field we cannot find 5 distinct points defined over F. 
Even for p = 101 the procedure will not choose distinct points in P 2 (Fioi) 
in about 1.5% of all cases. For p as large as 10007 the failure probability is 
neglectable. 

A method to get points which works for small finite fields is the follow- 
ing: Since points have codimension 2 in P 2 , we might use the Hilbert-Burch 
matrix, see |Ei95j . Theorem 20.15. We recall the Hilbert-Burch theorem in 
its most useful special case: 

Theorem 1.1. Let R be a local (or graded) regular noetherian ring. The 
minimal free resolution of R/I for a Cohen- Macaulay ideal I C R of codi- 
mension 2 has shape 

<- R/I <- R^ R n+1 <- R n <- 

where the map R <— R n+l is given by the n + 1 maximal minors of the matrix 
R n+1 R n . Conversely, given a matrix (cLij)i=o,...,n,j=i,...,n; whose maximal 
minors have no common factor, then the ideal generated by the maximal 
minors is Cohen- Macaulay of codimension 2 and the corresponding complex 
is exact. 



Note that (m , . . . , m n )(aij) 
panding the determinants 



= det 



for rrik 

/«oi • • • 
\a nl . . . 



= (— \) k det(eiy)i^fc follows by ex- 



of matrices with a repeated column with respect to the last column. The 
Hilbert-Burch Theorem is the basic reason why the Hilbert scheme of finite 



length subschemes on a smooth surface is smooth, see Fog68 
We continue with our Macaulay2-piogra.m 



19 : 
110 
ill 
il2 



M=random(S"{3:-4,l:-5},S"{3:-6}) ; — a Hilbert-Burch matrix; 
betti M 

IDelta=minors(3,M) ; — and its minors 
betti res IDelta 
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ol2 = total: 14 3 
0: 1 
1 
2 
3 

4: . 3 
5: . 1 3 

012 : BettiTally 

il3 : degree IDelta==delta 

013 = true 



STEP 3. Compute curves with double points in Delta: 

il4 : — the saturation of the ideal IDelta~2 contains 
— all equations which vanish double at Delta: 
J=saturate (IDelta~2) ; 

i!5 : betti J 



1 

ol5 = total: 1 10 
0: 1 
1 
2 
3 
4 
5 
6 
7 

8: . 1 
9: . 9 
ol5 : BettiTally 



As expected, there is precisely one curve with multiplicity 2 at every point 
of A. 



116 : IC=ideal (gens J*random (source gens J, S~{-d})) 

117 : degree IC 
ol7 = 9 
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Step 4. To certify that we have indeed obtained a 5-nodal curve C, it 
suffices to prove that the singular locus is reduced of degree 5 because only 
plane curves with at most nodes have a reduced singular locus. 

118 : singC=ideal jacobian IC + IC; 

119 : codim singC == 2 and degree singC == delta 
ol9 = true 

120 : — remove primary component at the irrelevant ideal : 

singCs=saturate(singC) ; 

121 : betti singCs, betti singC 



o21 = 



(total 

1 
2 
3 
4 
5 




1 
1 



1 

4, 



3 
1 



total 

1 
2 
3 
4 
5 
6 
7 



1 

1 3) 
1 



021 : Sequence 

i22 : codim (minors (2, jacobian singCs)+singCs) == 3 

022 = true 

We verify that C is geometrically irreducible from our information so far. 
Indeed, if the curve decomposes as C1UC2 with (degCq, degC^) = {a,b), say 
a < b, and a + b = 9, then the intersection points C\ fl C2 form a subset of 
A. So (a, b) — (1, 8) or (2, 7) is excluded because /a is generated in degree 6. 
The case (4, 5) is excluded since 20 > 18 = 5 and (3, 6) is excluded because 
A is not a complete intersection. 

Finally, we conclude from this computation that VJH W is unirational over Q. 
We first note that the computations in the finite prime field F may be viewed 
as the reduction mod p of the analogous computations for a curve defined 
over an open part of Spec Z. By semi-continuity, the curve over Q is 5-nodal 
as well, hence has geometric genus 10, which proves that we have a rational 
map A™ — » SDT10 defined over Q for a suitable n. More precisely, we have a 
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correspondence 



Hilb 5 (P 2 ) H= {(C, A) | C e FH°(F 2 ,l 2 A (d))* is 5-nodal } 

For g = 10 the left arrow gives a birational map of a component if' of H 
onto a dense open subset of Hilb 5 (P 2 ), because x(Zi(9)) = (") - 35 = 1, 
and /i°(2^(9)) = 1 holds for our specific point A G Hllb^P 2 ). 

The downward arrow factors through the universal C Vic d g — > 9Jl g , 
which has codimension at most (2 + l)(g + 2 — d) in the universal Picard 
variety over VJl g at every point. A fiber of H — > Wj over a point (C, L) with 
h°(C,L) = 3 is precisely the PGL(3) orbit. Thus to prove that if' ->■ 9Jti 
dominates, it remains to certify that the fiber of Wj — > Wlw over our specific 
point C has expected dimension p at the specific line bundle L = r]*0^2{\) e 
W$(C) where i] : C — > C C P 2 denotes the normalization map. 

By Riemann-Roch and adjunction h°(C,L) = r + 1 = 3 holds because 
h\C,L) = h (C,u c ® L- 1 ) = h°(¥ 2 ,l A (d - 4)) = 3. Moreover, the Petri 
map can be identified with 

H°(F 2 ,l A (d - 4)) ® H°(F 2 , 0(1)) -)■ H°(F 2 ,l A (d - 3)) 

which is an injection, since there is no linear relation among the three quintic 
generators of I a by the shape of the Hilbert-Burch matrix, compare also with 
the output line 08 and ol2 of the Macaulay2 program above. Thus WjiC) 
is actually smooth of dimension p in L as desired. 

It is easy to transform the code above into a function which chooses 
randomly 5-nodal curves of degree d provided the expected h°(F 2 ,I A (d)) = 
( d ~2 2 ) — 35 > 0. An implementation can be found in the Macaulay2 package 
RandomCurves. 

2 Searching 

Over a finite field F we might find points C G M(F) of a moduli space, 
provided M is dominated by a variety H of fairly low codimension in a uni- 
rational variety G. Indeed, if H is absolutely irreducible, then the proportion 
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of F-rational points is approximately 

I H(¥) 



G(¥) 



where q —\ F | and c is the codimension of H C G by the Weil formulas. If 
we can decide p ^ H fast enough computationally, then we might be able to 
find points in H(W q ) for small q by picking points at random in G(¥ q ) and 
testing p e H. I will illustrate this technique by searching for plane models 
of random curves of genus 11. Of course, to get just a random curve of genus 
11, it is much better to use Chang and Ran's unirational parameterization 
of OJtn via space curves, as indicated in the next section. 

This time we will use a bit more of the Macaulay2 syntax. In particular 
we will illustrate the use of method functions with options. 

i23: randomDistinctPlanePoints = method (TypicalValue=> Ideal); 

— create the ideal of k random points 
— via a Hilbert-Burch matrix 

i24: randomDistinctPlanePoints (ZZ ,PolynomialRing) := (k,S) -> ( 
if dim S =!= 3 then error "no polynomial ring in 3 variables' 

— numerical data for the Hilbert-Burch matrix 
n := ceiling((-3+sqrt(9.0+8*k))/2) ; 

eps := k - binomial (n+1 , 2) ; 

a := n+l-eps; 

b := n-2*eps; 

distinct := false; 

while not distinct do ( 

— the Hilbert-Burch matrix is B 
B := if b >= then random(S~a, S~{b : -1 , eps : -2}) 
else random(S"{a:0,-b:-l}, S~{eps:-2}); 

I := minors (rank source B, B) ; 

distinct = distinctPlanePoints I) ; 

return I) ; 

i25: distinctPlanePoints=method(TypicalValue=>Boolean) ; 
i26: distinctPlanePoints (Ideal) := I-> ( 

dim I==l and dim (I+minors(2, jacobian I))<=0) 
i27: distinctPlanePoints (List) : = L->( 

— for a List of ideals of points 

— check whether they have some point in common. 
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degree intersect L == sum(L, I->degree I)) 

The function randomDistinctPlanePoints will return an ideal of a set of 
distinct points which we use in our search below. 

The minimal degree of plane models of a general curve of genus g — 11 
is d — 10, and we expect that a general model will have 5 = ( d ^ 2 1 ) — 11 = 25 
ordinary double points. Since x(P 2 > ^a(IO)) = — 9 we expect that A with 
/i°(P 2 ,Zi(10)) = 1 from a codimension 10 subfamily of Hilb 25 (P 2 ). We can 
improve our odds if we look for models with 2 triple points. Since triple points 
occur in codimension 1 and p(ll, 2, 10) = 2 we expect that the family of plane 
curves of degree 10 with two triple points Pi,p 2 and 19 double points will 
dominate 5£R n . Since x(P 2 ,X| ® 0X^(10)) = ( 10 + 2 ) - 3 • 19 - 6 • 2 = -3, 
we expect a search for points in a subfamily of codimension 4. Thus the 
search function below will have an average running time of order 0(q A ) with 
respect to the number of elements of F = ¥ q . So this will be only feasible for 
very small finite fields. 

i28 : searchPlaneGenusllCurve=method(Options=>{Attempts=>inf inity}) 

— search a plane curve of degree 10, geometric genus 11 with 

— two triple points and 19 double points 

i29: searchPlaneGenusllCurve PolynomialRing := opt -> S -> ( 

11 := ideal(S_0,S_l) ; 

12 := ideal (S_l, S_2) ; 

IDelta := ideal 0_S; J := ideal 0_S; 
h := -3; attempt := 0; 

while h <= and attempt < opt. Attempts do ( 

IDelta = randomDistinctPlanePoints ( 19, S) ; 
while not distinctPlanePoints ({II, 12, IDelta}) do ( 
IDelta = randomDistinctPlanePoints ( 19, S) ) ; 
J = intersect (I1~3,I2~3, saturate (IDelta~2) ) ; 
h = (tally (degrees gens truncate (10, J) )_1)_{10}; 
attempt = attempt +1) ; 
print attempt; 

if attempt >= opt. Attempts then return null; 
f := ideal 0_S; 
gj := gens J; 

if h === 1 then f = ideal gJ_{0} 

else while f == do f=ideal (gJ*random (source gj, S~{ -10})); 
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singf := ideal singularLocus f; 
doublePoints := saturate (singf , 11*12); 
if degree doublePoints == 19 and ( 

f31 := ideal contract (S_2~7 , gens f ) ; 
dim singularLocus f31 == 1) 
and ( 

f32 := ideal contract (S_0~7 , gens f); 
dim singularLocus f32 == 1) 
then return f else 

searchPlaneGenusllCurve(S, Attempts => opt .Attempts-attempt)) ; 

After these preparations the commands below will return a desired curve 
within a few minutes. 

i30 : p=5;FF=ZZ/p — a finite ground field 

132 : S=FF[x_0. .x_2] 

133 : setRandomSeedC alpha") 

134 : C=searchPlaneGenusllCurve(S,Attempts=>2*p~4) 
18 

432 

— used 48.132 seconds 

o34 : Ideal of S 

3 Space Curves 

The proof of the unirationality of the moduli space Wl g for g = 11, 12 and 
13 by Sernesi and Chang-Ran is based on models of these curves in P 3 , sec 
jSer^lCRM] . 

Suppose C C P 3 is a Cohen-Macaulay curve with ideal sheaf Ic- The 
Hart shorne- Rao module 

M = Y,H\¥\Z c (n)), 

which has finite length and measures the deviation for C to be projectively 
normal, plays an important role in liaison theory of curves in P 3 . We briefly 
recall the basic facts. 

Let S = F[ Xq, . . . , X3] and Sc = S/Iq denote the homogeneous coordinate 
ring of P 3 and C C P 3 respectively. By the Auslander-Buchsbaum-Serre 
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formula, |Ei95] Theorem 19.9, Sc has projective dimension pd s Sc < 3 as an 
S'-module. Thus the minimal free resolution has shape 

<- S c <- S <- F l <- F 2 <- F 3 <- 

with free graded modules Fj = @S(— iY ij ■ 

The sheafified kernel Q = ker(F\ — > Oya) is always a vector bundle, and 

<- O c <- F s <- ©j<pp3(-j)^ <- G <- 

is a resolution by locally free sheaves. If C is arithmetically Cohen-Macaulay, 
then F3 = and Q splits into a direct sum of line bundles. In this case the 
ideal Ic is generated by the maximal minors of F\ <r- F2 by the Hilbert-Burch 
TheoremO In general we have M S J2 n& H 2 (F 3 , Q{n)) and Y. n & H 1 ^ 3 , Q{n)) = 
0. 

The importance of M for liaison comes about as follows. Suppose f,gE 
Ic are homogeneous forms of degree d and e without common factor. Let 
X = V(f,g) denote the corresponding complete intersection, and let C be 
the residual scheme defined by the homogeneous ideal Ic = (f,g) '■ Ic- 
The locally free resolutions of Oc and Oc are closely related: Applying 
£xt 2 {— , u v z) to the sequence 

-> Xc/x -^O x ^Oc^0 

gives 

^— £xt 2 (Z c /x, ^p 3 ) ^- ux <— ujc 0. 

From ujx — Ox(d + e — 4) we conclude £xt 2 (Ic/x, 0^z(—d — e)) = Cc*' and 
hence I^'pc — ^c(—d — e + 4). 
The mapping cone 

0^O c ^<^ ©OH)^. 

Ox O - O(-d) © 0(-e) 0(-d - e) 

dualized with "Hom(— , 0(—d — e)) yields the locally free resolution 

-> 0{j - d - efv -> - e) © O(-e) © O(-d) -> -> G>c -> 
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In particular one has 



M c , = H 1 ^, x o(n)) = H 1 ^ 3 , Q*{n-d- e)) 

= H 2 (W 3 , g(d + e-4- n))* ^ Hom F (M c , F)(4 - d - e) 

ngZ 

Thus curves, which are related via an even number of liaison steps, have the 
same Hartshorne-Rao module up to a twist. Rao's famous result [Rao78j 
says that the even liaison classes are in bijection with finite length graded 
5- modules up to twist. 

Reversing the role of C and C, we conclude that the ideal sheaf of C has 
a locally free resolution 

<- l c <- T © Ci <- £ 2 <- 

where L\ = ®iO(—ci) and £ 2 = Q)kO(—d k ) are direct sums of line bundles, 
while J 7 is a locally free sheaf without line bundle summands. Note that T 
has no i7 2 -cohomology, and its H ^-cohomology 

is the Hartshorne-Rao module of C. The map £2 £1 coming from a liaison 
construction might be non-minimal, in which case one can cancel summands 
of £1 against some summands of C<i- Below we will work with complexes 
which arise after such cancellation. 

Since there are no line bundle summands and no H 2 , T is determined by 
its if ^-cohomology: Consider a minimal presentation of M 

<- M 4r- ® l S(-a l ) <- ®jS{-bj) <- N <- 

and the kernel N. Then T = N is the associated coherent sheaf of N. 

The main difficulty in constructing space curves lies in the construction 
of M. Given M we can find the curve as the cokernel of an homomorphism 
<f G Hom(£ 2 ,^ r © £1), which we may regard as a vector bundle version of 
a Hilbert-Burch matrix. Frequently in interesting examples we will have 
£1 = 0. 

We apply this approach for the case g = 12 and d = 13. 



12 



To get an idea about the twists in the resolutions of M and C, we use 
the Hilbert numerators: Suppose the minimal finite free resolution of M is 



<- M 



F <- F 1 <- . . . <- F 4 <- 



with 

then the Hilbert series 

H M (t) = ^ dim MJ 



Assuming the open condition that C has maximal rank, i.e. 

H {F 3 ,O(n)) ^H°{C,L n ) 

is of maximal rank for all n then Hj^it) = 5t 2 + 8t 3 + 6t 4 and the Hilbert 
numerator is 

H M (t)(l - t) 4 = ht 2 - 12t 3 + 4t 4 + 4t 5 + 9t 6 - lQt 10 + 6t n 

If M has a natural resolution, which means that for each j at most one 
is non-zero, then M has a Betti table 



total 
2 
3 
4 



12 3 4 

5 12 17 16 6 

5 12 4 . . 

. . 4 . . 

. . 9 16 6 



in Macaulay2 notation. Having natural resolution is another open condition. 
If Sc = S/Ic has natural resolution as well, then its Betti table will be 



total 

1 
2 
3 
4 
5 



12 3 

1 11 16 6 
1 . 



2 

9 16 6 



13 



Comparing these tables, we conclude that C\ = 0, rank J 7 = 12 — 5 = 7 and 
£-2 = 0(— 4) 4 © 0{— 5) 2 with rank£ 2 = 6 will hold for an open set of curves. 
Of course, we have to prove that the set is non-empty and contains smooth 
curves. 

To construct a desired M, we start with the submatrix tp defining S 12 (— 3)<— S 4 (— 4) 
in F\ -f— F2 which we choose randomly. Since 12 ( 3 ^ 1 ) — 4( 3 ^ 2 ) = 8 > 5, the 
kernel ker('0* : S' 12 (3) — > S 4 (4)) has at least 8 linearly independent elements 
of degree —2, and we choose 5 random linear combinations of these elements 
to get the transpose of the presentation matrix F <— F\ of M. 

135 : FF=ZZ/10007 

136 : S=FF[x_0. .x_3] 

137 : psi=random(S~{12:-3},S~{4:-4}) — the submatrix 

138 : betti (syzpsit=syz transpose psi) 

139 : M=coker transpose (syzpsit*random(source syzpsit , S~{5 : 2}) ) ; 

140 : F= res M — free resolution of the desired Hartshorne-Rao module 

141 : betti F 

142 : L2=S~{4:-4,2:-5} 

143 : phi= F.dd_2*random(F_2,L2) ; 

144 : betti (syzphit=syz transpose phi) 

145 : IC=ideal mingens ideal ( transpose syzphit_{5}*F . dd_2) ; 

146 : betti res IC — free resolution of C 

147 : codim IC, degree IC, genus IC 
o47 = (2, 13, 12) 

o47 : Sequence 

Next we check smoothness. Since CcP 3 has codimension 2 and is locally 
determinantal, it is unmixed, and we can apply the jacobian criterion. 

148 : singC=IC+minors (2, jacobian IC) ; 

149 : codim singC==4 
o49 = true 

Thus the curve C is smooth. We conclude that the Hilbert scheme Hilbd iS (P 3 ) 
of curves of degree d — 13 and g — 12 has a unirational component H defined 
over Q, whose general element is smooth. 

As a corollary we get the unirationality of 97ti2 from this if we verify 
that the fiber of the rational map H --■» SDT12 has the expected dimension 
dimPGL(4) + p(g,r,d) = 15 + 4 at the given point. Now this holds if 
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the Petri map H (u) C (—l)) ®H (O(l)) ->■ H°(u c ) is injective where u c = 
£xt 2 {Oci^v z ) — Sxt l (Zc, Ofz{— 4)) denotes the canonical bundle on C. 

i50 : betti Ext"l(IC,S~{-4}) 
1 

o50 = total: 6 12 
-1: 2 . 
0: 4 12 

shows that there are no linear relation among the two generators of 
r,(w c ) = ^Ti^PVcW) = Extl(I c ,S(-A)) 

in degree —1. Thus H dominates and Wl±2 is unirational. As a corollary of 
our construction we obtain that the Hurwitz scheme Hk,g °f fc-gonal curves 
of genus g is unirational for (k,g) = (9, 12). Indeed wc(- 1) is a line bundle 
of degree fc = 22 - 13 = 9. 

The case d = 13, g = 12 is actually not used in Sernesi's proof of 
the unirationality of Wlw Chang, Ran and Sernesi use the cases (d, g) = 
(10, 11), (12, 12) and (13,13), all of which can be treated similarly as the 
case above. I took the case (d, g) = (13,12) because it illustrates the dif- 
ficulty in constructing M well, and because being different, it implies some 
minor new results. Computer algebra simplifies the cumbersome proof that 
the approach leads to smooth curves in all of theses cases. For more details 
and an implementation we refer to the Macaulay2 package RandomCurves. 

There are 65 values for (d, g) such that there are possibly non-degenerate 
maximal rank curves with natural resolution such that the Hartshorne-Rao 
module has diameter < 3 and natural resolution as well. For 60 of these 
values one can establish the existence of a unirational component in the 
Hilbert scheme by the methods above. 

To prove the unirationality of VJl g for g > 14 via space curves, leads to 
Hartshorne-Rao modules which have diameter > 4, i.e., modules, which are 
nonzero in at least 4 different degrees. The construction of such modules is 
substantially more difficult. 

4 Verra's proof of the unrationality of DJlu 

Verra's idea [Ve05j is of beautiful simplicity. Consider a general curve of 
genus g = 14 and a pencil \D\ of minimal degree which is degD = 8. The 
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Serre dual linear system \K — D\ embeds C into P 6 as curve of degree 18 
with expected syzygies 

1 2 3 4 5 
total: 1 13 45 56 25 2 

0: 1 

1: . 5 . . . . 
2: . 8 45 56 25 . 
3: 2 

In particular, C C P 6 lies on 5 quadrics. The intersection of the quadrics 
should have codimension 5 and degree 2 5 = 32, thus should equal CUC" 
where the residual curve C has smaller degree 32 — 18 = 14 than C and 
hence also smaller genus, g{C') = 8 as it turns out. 

(Indeed the dualizing sheaf of the complete intersection is ucuc — Ccuc*'(5' 
2-7) = 0cuc(3), and we obtain arithmetic genus p a {CUC') = 2 5 -3/2 + 1 = 
49. Assuming only nodes as intersection points we have ucuc ® @c — 
u c {C n C) and get deg(C n C) = 3 ■ 18 - 26 = 28 intersection points. 
The formula p a (C U C) = p a {C) + p a {C) + deg(C n C) - 1 finally gives 
Pa{C) = 49 — 14 — 28 + 1 = 8.) There is no reason to expect anything else 
than, that C is a smooth curve. 

By a famous result by Mukai |Mu92j, every general canonical curve of 
genus 8 is obtained as transversal intersection 

C = G(2, 6) n P 7 C P 14 

of the Grassmannian G(2, 6) C P 14 in its Pliicker embedding. If we choose 
8 general points on G(2, 6) and P 7 as their span, then we get a genus 8 
curve together with 8 points. We group these 8 points into two divisors 
D\ — pi + . . . +p4 and D 2 = p$ + • • • +ps of degree 4. Then | Kc> + D\ — D 2 \ 
is a general linear system of degree 14 on C, and re-embedding C ■=->■ P 6 
leads to a curve with expected syzygies 
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The Macaulay2 code is now straightforward. First we construct C in its 
canonical embedding. 
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i51 : randomCurveOf Genus8With8Points = R ->( 

— Input: R a polynomial ring in 8 variables, 

— Output: a pair of an ideal of a canonical curve C 

together with a list of ideals of 8 points 
— Method: Mukai's structure theorem on genus 8 curves. 

Note that the curves have general Clifford index. 
FF:=coeff icientRing R; 
p:=symbol p; 
— coordinate ring of the Plucker space : 
P:=FF [flatten apply (6, j->apply(j , i->p_(i, j)))] ; 
skewMatrix:=matrix apply(6,i->apply(6, j->if i<j then 
p_(i,j) else if i>j then -p_(j,i) else 0_P)); 

— ideal of the Grassmannian G(2,6): 
IGrass : =pf af f ians (4, skewMatrix) ; 

points : =apply (8 , k->exteriorPower (2 , random (P "2 , P~6) ) ) ; 
ideals : =apply(points ,pt->ideal ( vars P*(syz pt**P~{-l}) ) ) ; 

— linear span of the points: 

L:= ideal (gens intersect ideals) _{0. .6}; 

phi:=vars P°/ L; — coordinates as function on the span 

— actually the last 8 coordinates represent a basis 
phi2:= matrix-ftoList (7 : 0_R)} I vars R; 

— matrix for map from R to P/IC 

IC:=ideal (gens IGrass°/„L) ; — the ideal of C on the span 

— obtained as the reduction of the Grassmann equation mod L 
IC2:=ideal mingens substitute (IC, phi2) ; 

idealsOf Pts : =apply (ideals , Ipt-> 

ideal mingens ideal sub (gens Ipt%L,phi2) ) ; 
(IC2, idealsOf Pts)) 

Building upon Mukai's result, we can construct the desired curve C: 

i52 : randomNormalCurveOf Genus8AndDegreel4 = S -> ( 

— Input : S coordinate ring of P~6 

— Output: ideal of a curve in P~6 
x:=symbol x; 

FF:=coeff icientRing S 
R:=FF[x_0. .x_7] ; 

(I, points) : =randomCurveOf Genus8With8Points (R) ; 

Dl : =intersect apply (4, i->points_i) ; — divisors of degree 4 
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D2:=intersect apply (4, i->points_ (4+i) ) ; 

— compute the complete linear system I K+D1-D2 I , note K=H1 
Hl:=gens Dl*random(source gens Dl , R~{-1>) ; 

El :=(I+ideal H1):D1; — the residual divisor 
L:=mingens ideal(gens intersect (El ,D2)°/„I) ; 

— the complete linear system 

— note: all generators of the intersection have degree 2. 
RI:=R/I; — coordinate ring of C in P~7 

phi : =map (RI , S , substitute (L , RI) ) ; 
ideal mingens ker phi) 
i53 : FF=ZZ/10007;S=FF[x_0. .x_6] ; 

155 : I=randomNormalCurveOf Genus8AndDegreel4(S) ; 

156 : betti res I 

Finally, we get curves of genus 14 with 

i57 : randomCurveOf Genusl4=method(TypicalValue=>Ideal , 
Options=>{Certif y=>f alse}) 

— The default value of the option Certify is false, because 

— certifying smoothness is expensive 

158 : randomCurveOf Genusl4(PolynomialRing) :=opt ->( S-> ( 

— Input : S PolynomialRing in 7 variables 

— Output : ideal of a curve of genus 14 

— Method: Verra's proof of the unirationality of M_14 
IC :=randomNormalCurve0fGenus8AndDegreel4(S) ; 

— Choose a complete intersection: 

CI:=ideal (gens IC'*random (source gens IC , S~{5:-2})) ; 
IC:=CI:IC; — the desired residual curve 
if not opt . Certified then return IC; 

if not (degree IC ==18 and codim IC == 5 and genus IC ==14) 

then return nil; 
someMinors :=minors(5, jacobian CI); 
singCI : =CI+someMinors ; 

if not (degree singCI==28 and codim singCI==6) 

then return nulll; 
someMoreMinors:=minors(5, jacobian (gens IC)_{0. .3,5}) ; 
singC : =singCI+someMoreMinors ; 

if codim singC == 7 then return IC else return nil)) 

159 : time betti ( J=randomCurveOf Genusl4(S) ) 
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160 : time betti(J=randomCurveOfGenusl4(S,Certif ied=>true) ) 

161 : betti res J 

To deduce from these computation the unirationality of Wtu, we have to 
prove again that the Petri map is injective. By the Betti numbers of the 
free resolution of Sc we see that there is is no linear relation among the two 
generators of u>c(—l)- Thus the family is dominant because the conditions 
that C has expected Betti numbers and that C is smooth are open. 

Finally, we remark that, comparing the syzygies of Oc, OcuC and 0& 
via liaison as outlined for space curves in Section [3j we see that the Petri 
map of (C, Oc(l)) is injective, if and only if la is generated by quadrics. 
Indeed, suppose the Betti table of C is 



where we assume that we need x cubic generators of the ideal of la- Then 
the Betti table of the mapping cone with the Koszul complex resolving Ocuc 
is 



1 2 3 4 5 6 

-1: . 1 

0: 1 . 5 . . . . 

1: . 7 x 10 . . . 

2: . x 35 56 45 8 . 

3: 5 . 

4: 1 



Minimalizing the dual complex, leads to following Betti table of C: 
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8 45 56 25 x 
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5 Minimal Resolution Conjectures and Koszul 
Divisors 



A graded ^-module M is said to satisfy the minimal resolution conjecture 
(MRC) (or is said to have expected syzygies) if the Betti numbers of the 
minimal free resolution 

* M F - Fi - ...F c -> 

where Fi = (BjS(— jY ij satisfy the following: For each internal degree j at 
most one of the numbers (3ij ^ 0. In other words, M satisfies MRC if one 
can nearly read off the Betti table from the Hilbert numerator of the Hilbert 
series Hm- 

We say that a module M has a pure resolution, if for each homological 
degree % at most one /3y ^ 0. Betti tables of Cohen- Macaulay modules with 
pure resolutions play a special role since they span the extremal rays in the 
Boij-Soderberg cone of all Betti tables, see |ES09[lBS08j . Moreover, if M has 
a pure resolution, then M satisfies the MRC, see [ES09] . 

We say that the MRC holds generically on a component H of a Hilbert 
scheme Hilb p ( t )(P n ) if it is satisfied for the coordinate ring Sx for X G U C 
H C Hilb p ( t )(P n ) in a dense open subset U of H. Note that the Hilbert 
function is constant on an open set U' of H, so this makes sense. Since Betti 
numbers are upper semi- continuous in flat families with constant Hilbert 
function, there is a smallest possible Betti table for X 6 U', and we can ask 
what this table is. If the MRC is satisfied generically on H, then in a sense, 
the question, what is the generic Betti table, has a trivial answer. 

The MRC has been established in various cases, the most famous one is 
the following. 

Theorem 5.1 (Voisin, 2005 jVo02] . Generic Green's Conjecture |Gr84] ). A 
general canonical curve of genus g > 3 over a ground field of characteristic 
satisfies the MRC. 

In more concrete terms, this conjecture says for example, that a general 
canonical curve of genus g = 15 has the Betti table 

1 2 3 4 5 6 7 8 9 10 11 12 13 
total: 1 78 560 2002 4368 6006 4576 4576 6006 4368 2002 560 78 1 
0: 1 
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1: . 78 560 2002 4368 6006 4576 

2: 4576 6006 4368 2002 560 78 . 

3: 1 

Note that the symmetry of the table comes from the fact that the homo- 
geneous coordinate ring of a canonical curve is Gorenstein. Before Voisin's 
result, g = 15 was for a long time the bound of how far a computer algebra 
verification of the generic Green Conjecture was feasible, the main problem 
being the memory requirements in the syzygy computation. 
The full Green conjecture is the following 

Conjecture 5.1 (Green's Conjecture |Gr84j ). Let C be a smooth projective 
curve of genus g > 3 over a field of characteristic zero. The canonical ring 
Rc = s ^ n H°{C 1 ujf; n ) of C as an S = SymH°(C,coc) -module has vanishing 
Betti numbers f3 iii+2 = if and only if i < Cliff (C) . 

Here the Clifford index is defined as 

Cliff(C) = min{degL-2(/i°(L)-l) | h°{L),h\L) > 2}. 

This conjecture generalizes the Noether-Petri-Babbage theorem which are 
the cases i = and % = 1 respectively. Its known for i = 2 in full generality, 
see |Vo88t ISch91j . The direction Cliff (C) < % =>- /3j,j+2 ^ was established 
by Green and Lazarsfeld in the appendix to |Gr84j 

Green's conjecture is known to be wrong over fields of finite characteristic. 
The first cases are g = 7 in characteristic 2 and genus g = 9 in characteristic 
3, see |Sch861ISch03] . 

Frequently, Green's Conjecture is formulated with the property N p of 
Green-Lazarsfeld |GL85j : A subscheme X C P r satisfies the property N p if 
Pij = for all j > % + 2 and i < p, in other words if the first p steps in the free 
resolution are as simple as possible. An equivalent formulation of Green's 
Conjecture is that a canonical curve satisfies property N p iff p < Cliff (C). 

If the coordinate ring Sx for general X G U' C H C Hilb p ( t )(P r ) is pure, 
then 

Z = {X G U' | S x does not satisfy MRC} 

is a divisor, a so-called Koszul divisor, because in principle it can be computed 
with the Koszul complex and Koszul cohomology groups. 

For odd genus g = 2k — 1 Hirschowitz and Ramanan computed the class 
of the corresponding divisor in 97t r 
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Theorem 5.2 (Hirschowitz-Ramanan [HR98J). If the generic Green conjec- 
ture holds for odd genus g = 2 k — 1, then the Koszul divisor 

Z = {C eTi g \ R c does not satisfy MRC } = (k - 

where 9Jt* fc = {C G VJt g \3gl on C} is the Brill-Noether divisor. 

The proof is based on a divisor class computation of Z and DJl g k on a 
partial compactification of 9Jt 9 and the fact that Z — (k — l)9Jt^ fc is effective. 
The coefficient A; — 1 is explained by the fact that a curve with a g\ has 
0k-i,k = Pk-2,k > k — 1 and equality holds for a general fc-gonal curve of 
genus g = 2k — 1 as a consequence of the Hirschowitz-Ramanan Theorem. 

Based on the Hirschowitz-Ramanan and Voisin's Theorem, Aprodu and 
Farkas |AF09j established that Green's Conjecture holds for smooth curves 
on arbitrary K3 surfaces. 

One can hope that there are various further interesting Koszul divisors. 
To get divisors on Wl g we consider the case that g = (r + l)s, d = r(s + 1) 
and ask whether a normal curve of degree d, genus g and speciality h l (L) = s 
can have a pure resolution. This is a purely numerical condition on r and s. 

Conjecture 5.2 (Farkas, |Fa06j ). Let p > and s > 2 be integers. Set 
r = (s + l)(p + 2) — 2, g = (r + l)s and d = r(s + 1). A general smooth 
normal curve C C P r of genus g, degree d and speciality h x {OciX)) = s ^ as 
a pure resolution, equivalently it satisfies N p . 

If the conjecture is true, then the divisor 

Zp, s = {C £ I which does not satisfy property N p } 

gives counterexamples to the slope conjecture, see [AF08] . For small cases 
Farkas verified the conjecture computationally in the spirit of the computa- 
tions below. 

Turning to non-special curves, we have 

Conjecture 5.3 (Farkas, [FLIP] ). A general smooth normal curve C C P r 
of odd genus g = 2p + 3 and degree d = 2g has a pure resolution, equivalently 
satisfies N p . A general smooth normal curve of even genus g = 2p + 6 and 
degree 2g — 2 has a pure resolution, equivalently satisfies N p . 
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For the rest of this section we return to the task to verify these conjectures 
computationally in a few cases. The case of genus g = 8 and degree d — 14 
has been established along the proof of Verra's Theorem in Section HI We 
treat the case of genus g = 7 and degree d = 14. This case was settled in 
[Fa09j using reducible curves. Here we use smooth curves. 

Consider for g = 7 a random plane model of degree d = 7 with 5 = 8 
nodes. 

i62: FF=ZZ/10007 

i63: R=FF[x_0. .x_2] 

i64: g=7 

i65: delta=binomial(6,2)-7 

i66 : J=randomDistinctPlanePoints(delta,R) 

i67: betti res J 

i68: betti (J2=saturate(J~2) ) 

i69: C=ideal (gens J2*random (source gens J2,R~"C _ 7})) 

To find a general divisor of degree 2g = 14 on C we note that over a 
large finite field geometrically irreducible varieties have always a lot of F- 
rational points. We can frequently find points as one of the ideal theoretic 
components of the intersection with a random complementary dimensional 
linear subspace. 

i70: decompose (C+ideal random (1,R)) 

i71 : apply (decompose (C+ideal random(l ,R) ) , c->degree c) 
i72: time tally apply ( 1000, i-> 

apply (decompose (C+ideal random(l ,R) ) , c->degree c)) 
— used 9.62383 seconds 
o72 = Tally-C-Cl, 1, 1, 1, 1, 2} => 5} 

{1, 1, 1, 1, 3} => 18 

{1, 1, 1, 2, 2} => 21 

{1, 1, 1, 4} => 47 

{1, 1, 2, 3} => 84 

{1, 1, 5} => 91 

{1, 2, 2, 2} => 19 

{1, 2, 4} => 120 

{1, 3, 3} => 58 

{1, 6} => 167 
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{2, 2, 3} => 40 
{2, 5} => 94 
{3, 4} => 77 
{7} => 159 
o72 : Tally 

Thus the following function will provide points. 

i73 : randomFFRationalPoint=method() 

i74: randomFFRationalPoint (Ideal) :=I->( 

— Input: I ideal of a projective variety X 
— Output: ideal of a FF-rational point of X 
R:=ring I; 

if char R == then error "expected a finite ground field"; 
if not class R === PolynomialRing then 

error "expected an ideal in a polynomial ring"; 
n:=dim 1-1; 

if n==0 then error "expected a positive dimensional scheme"; 
trial :=1; 
while ( 

H:=ideal random(R~l,R~{n:-l}) ; 
pts : =decompose (I+H) ; 
ptsl : =select (pts ,pt-> degree pt==l and dim pt ==1); 
#ptsl<l ) do (trial=trial+l) ; 
ptsl_0) 

i75: randomFFRationalPoint (C) 

This allows to get a random effective divisor of degree 14 on C with all 
points F-rational, for which we verify N 2 : 

i76 : points=apply(2*g, i->randomFFRationalPoint (C) ) 

i77: D=intersect points — effective divisor of degree 14 

i78: degree D 

o78 = 14 

i79: D J=intersect (D , J) 

i80 : degree D J==degree D + degree J 

o80 = true 

i81: betti DJ 
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i82: H=ideal (gens DJ*random( source gens DJ, R~{-6}))+C 
i83: E=(H:J2):D — the residual divisor 
i84: degree E +degree D + 2*degree J == 6*7 
o84 = true 

i85: L=mingens ideal ((gens truncate (6, intersect (E, J))) °/„C) 

— the complete linear series |D| 
i86: RC=R/C 
i87: S=FF[y_0. .y_7] 
i88: phi=map(RC,S,sub(L,RC)) 

i89: I=ideal mingens ker phi; — C re-embedded 
i90: (dim I, degree I, genus I) == (2,14,7) 
i91: time betti res I 

1 2 3 4 5 6 
o91 = total: 1 14 28 56 70 36 7 

0: 1 

1: . 14 28 . . . . 

2: . . . 56 70 36 7 

The computation above is the reduction mod p of the computation for a 
curve together with a divisor over some open part of Spec Ok of a number 
field K. We can bound the degree [K : Q] < 7 14 of the number field, a bound 
on its discriminant seems out of reach. 

Hence as before we can conclude that Farkas' syzygy conjecture holds 
for curves of genus 7 and degree 14 over fields of characteristic zero. Using 
Mukai's description of curves of genus 7 |Mu92j, we could find an example 
defined over Q. I chose to present the example above because I wanted to 
illustrate the trick to get points in Wl g , n from points in Wl g over finite fields. 

It is not difficult to use similar constructions to verify Farkas' conjecture 
for other small g. 

It is known that the MRC is not satisfied in a number of important cases, 
notably for non-special normal curves of genus g > 4 and large degree [Gr84j . 
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